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Interfaces determined by capillarity and gravity in a 
two-dimensional porous medium 

M. Calle* C. M. Cuestaj J. J. L. Velazquez^ 


Abstract 

We consider a two-dimensional model of a porous medium where circular grains 
are uniformly distributed in a squared container. We assume that such medium is 
partially filled with water and that the stationary interface separating the water phase 
from the air phase is described by the balance of capillarity and gravity. Taking the 
unity as the average distance between grains, we identify four asymptotic regimes that 
depend on the Bond number and the size of the container. We analyse, in probabilistic 
terms, the possible global interfaces that can form in each of these regimes. In 
summary, we show that in the regimes where gravity dominates the probability of 
configurations of grains allowing solutions close to the horizontal solution is close to 
one. Moreover, in such regimes where the size of the container is sufficiently large 
we can describe deviations from the horizontal in probabilistic terms. On the other 
hand, when capillarity dominates while the size of the container is sufficiently large, 
we find that the probability of finding interfaces close to the graph of a given smooth 
curve without self-intersections is close to one. 
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1 Introduction 

In this paper, we present a two-dimensional model of a porous medium that consists of a 
squared container partially hlled with a liquid and where the pore throats are formed by 
hxed circular grains that are uniformly distributed. We analyse in probabilistic terms the 
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structure of the stationary interfaces that separate the liquid from air. We consider several 
parameter regimes where, essentially, either gravity or capillarity dominates. 

In this model the form of the interfaces, both at the pore scale and at the macro scale, 
is determined by the balance between gravity and capillary forces. For a given volume of 
liquid, these interfaces are the union of single interfaces (see 0) that meet the solid matrix 
at given contact angles following Young’s law (see e.g. [I])- 

The quasi-stationary system seems to evolve by means of a sequence of local minimiz- 
ers of the energy. The energy for steady states is the sum of surface energy and potential 
gravitational energy. The contact angle conditions can be encoded in the different inter¬ 
facial energies. One of the most relevant features in these flows is the existence of a large 
number of equilibria that are local minimizers of the total energy of the system. In quasi- 
stationary situations, if the volume of the liquid phase changes slowly, transitions occur 
between different minimizers. Understanding the evolution of such a system would require 
the description of the global motion of the interface that explains the jumps between these 
multiple local energy minimizers. A hrst step towards this goal consists in understanding 
the structure of the minimizers. 

The aim of this paper is thus to describe the structure of equilibrium interfaces for 
different parameter regimes in the model described above. There are several characteristic 
lengths in the problem, namely, the distance between grains, their size, as well as the 
macroscopic size of the system. On the other hand, there are other length scales (the 
capillary length and the inverse of the Bond number) that, as we shall see, give the distances 
where surface tension and gravity balance (typically, surface tension dominates for small 
distances and gravity for large ones). 

This study is motivated by the experiments presented in Furuberg-Maloy-Feder m 
and [12] • The authors reproduced Haines jumps (see [S]) in a slow drainage process in 
a container consisting of two plates separated by thin cylinders that act as obstacles or 
‘grains’. We recall that Haines jumps are abrupt changes of the pressure in the liquid phase 
due to sudden changes in the geometry of the interface. From the mathematical point of 
view these jumps seem to correspond to a fast transition from one equilibrium state to a 
different one in the capillarity equations. 

We mention that there are several mathematical studies related to the study of Mul¬ 
tiphase Flow in Porous Media. However, much of the mathematical analysis has focused 
on the qualitative study of semi-empirical macroscopic descriptions of these flows (see e.g. 
0)- There are also rigorous results concerning the derivation of such macroscopic laws 
from the flows taking place at the pore scale. Not surprisingly, such attends are restricted 
to single-phase flow (see the chapter of [13], [2D]) or simplihed models mimicking essential 
features of multi-phase flow (e.g. i, ng, na, Hi). Although our model falls into the 
second category, we believe that one can analyse simple, but not phenomenological, models 
in probabilistic terms to get a better understanding of the effect that processes at the pore 
scale have in the whole system. In this regard, our work is also related to for example [1] 
(and the references therein), where the propagation of fronts through a random held of 
obstacles is analysed in relation to hysteresis. 

In the next section we set up the mathematical model and discern the parameter regimes 
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under consideration. We outline the contents of the paper at the end of this Introduction. 


1.1 The mathematical model and main results 

We consider a two dimensional rectangular container of length L that contains N uniformly 
distributed circular grains. If the average distance between grains is d > 0 (understood as 
the average of the distances of every center to its nearest neighbour), then N oc /d?. We 
then let > 0 represent the density of grains, i.e. v satishes N = jd?. We assume, for 
simplicity, that all grains are of equal size and we denote by i? < d/2 their radius. 

The container is partially hlled with a liquid so that an interface forms separating the 
liquid from the air. We assume that such an interface splits the domain in two regions, i.e. 
it does not intersect itself and there are no isolated bubbles of either of the fluid phases. 
Let us denote by T the curve that joins grains and forms the interface, that is the union 
of capillary (-pressure) curves 7 joining every two grains that meet the interface and that 
we call elemental components of T. A sketch of this setting is shown in Figured] 



L 


Figure 1 : Schematic picture showing the problem set-up and the three relevant length 
scales. 

Let us £x the canonical basis {ei,e2} such that ei = ( 1 , 0 ) and 62 = ( 0 , 1 ), then 62 
points in the opposite direction to the gravitational held. We shall denote the coordinates 
of a point x G by x = (xi, 0:2) in this basis. 

Every elemental component of F, 7, satishes the capillary equation (e.g. m) 

aH = pgx 2 -p, 

that must be complemented with Young’s condition at the grains and/or the walls of the 
container. We shall assume that the contact angle is the same at all solid surfaces. Here 
H denotes the signed curvature of each 7 that forms F, p is the (constant) density of the 
liquid, g is the gravity constant and a denotes surface tension. Here p is the pressure 
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associated to the volume fraction of the domain (including grains), vi G (0,1), occupied 
by the liquid. Then 

p = pgX with X = viL . 

We take as unit of length the average distance between the grains of the system, d. Up¬ 
scaling the length variables with d, we obtain that T satishes the non-dimensional equation 

H = B{x2-X). (1) 


subject to a contact angle boundary condition that we specify below. We obtain the 
non-dimensional number B, which is the Bond number 


B 
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and measures the balance between gravitational forces and surface tension in this units of 
length. 

For simplicity, we use the same notation for the non-dimensional variables. Also we 
let L and R be the non-dimensional size of the domain (i.e. L/d) and grain radius (i.e. 
R/d), respectively. This, in particular, means that, with this notation, A = viL in ([1]). 
We observe that Vi depends on the conhguration of grains, so if we let Vq G [0,1] be the 
volume fraction of the domain occupied by the liquid, the volume occupied by the liquid 
is and then 

ViL^ VqL^ + miR^ (3) 

is the volume below the interface T, if n is the number of grains below T. 

We assume thereafter that Vq is hxed, and we take the horizontal line = ViL{= A) 
as a reference height for each conhguration of grains. Observe that this is a solution of ([1]) 
with a = 0 at the walls. 

We also observe that with the new notation, where L and R are non-dimensional 
lengths, then the number of grains in the container is Henceforth, with this 

scaling we can think of the equivalent problem where the liquid and air occupy some space 
in a square domain [0,L]^ with L ^ 1, where vLF' grains are uniformly distributed with 
an average distance 1 between them. 

We indicate the position of a center by the variable ^ ^ ^ind write its 

boundary in the form dB^{^). We take n as the unit normal vector to T pointing in the 
direction of air, and t the tangent vector to T such that {t, n} has the same orientation as 
the canonical basis (i.e. that det(t,n) = 1 ). 

Let us denote by N the outer normal vector to dB^l^) at a point P where T and dB^l^) 
intersect (i.e. P is a contact point). Let us denote by tc the value of t at P. Then, the 
Young condition is dehned as follows: if N ■ tc> 0, 


N A tc = sin(a)e3 (4) 

where a G [— f , f] is the contact angle and where 63 = ei A 62. If N • tc < 0 , then with 
the same choice of ei, we set N A tc = sin(a — 7r)e3. 
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We shall distinguish different parameter regimes. The following observations are then in 
order. We recall that B is the Bond number with respect to the scale given by the distance 
between grains. In the case 5 S> 1 the gravitational effects are the dominant ones. In this 
limit the liquid-gas interfaces tend to be horizontal lines, corrected by boundary layers of 
height In the limit i? <C 1, the interfaces are dominated by surface tension and this 

allows vertical variations in the connection between grains (of the order of the distance 
between the grains). Thus, unlike in the case i? S> 1, if i? <C 1 the connections can take 
place in any direction, not just horizontally. 

The above regimes take into account the size of B only, but we shall relate it to L as 
well. In particular, the limit L ^ oo gives an ever increasing number of grains and can in 
principle give rise to ’anomalous’ situations (not sketched above) that we analyse in detail 
in Section [3] and that we summarise below. 

We distinguish several cases depending on the relative size of B and L. In the cases 
where i? 3> 1, we actually compare L to min{A/B, which relates the height of the 

boundary layers to R: 

(i) The regime 1 < ^ < mm{y/B,R where 5 —)■ oo and i? —)■ 0 (Section 13.11 

Theorem 13.3p : In this case gravity dominates and the size of the container is not 
sufficiently large to allow signihcant capillarity effects. We hnd that as i? —)■ oo there 
are configurations with probability close to 1 for which the solution interfaces tend to 
the horizontal solution. One can in fact distinguish the cases L = 0(1) and L —)■ cx), 
with either R 1/\fB or i? ^ 1/ \fB, with the same result. 

(ii) The regime L = 0(min{\/5, i?“^}) and L —)■ oo (Section 13.21 Theorem [3.511 : 

In this case we find that the probability of finding solutions close to the horizontal 
line is close to one. On the other hand we also hnd that the probability of hnding 
conhgurations that have a given number of elemental components is close to zero, 
but bounded from above by a Poisson distribution. Again, the distinction between 
the cases l/\fB ^ R and l/y/B i? is not signihcant. 

(iii) The regime L min{\/5, and L —)■ oo: 

(a) The case B ^ oo (and i? —)■ 0) (Section 13.31 Theorem I3.10p . This is a 
regime where gravity in principle dominates, we hnd an interval of L in terms 
of B where with probability close to one solutions stay in a narrow strip around 
the horizontal reference line. For L above this range there is no guarantee that 
solutions are necessarily close to the horizontal reference line. If l/\fB 3> R, 
the interval is y/B L -C \fB\og{B), and if 1/y/B <C R with i? —)■ 0, the 
result holds replacing 1 / y/B by R. 

(b) The case B = 0{1). In this case, capillary and gravity forces are comparable 
and it is not possible to derive general properties about the structure of the 
solution interfaces using perturbative methods. Therefore this regime is not 
included in our analysis. 
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(c) The case 5 —)■ 0 (Section 13.41 Theorem 13.121) . We will see that in this case 
the only meaningful limit is that of L = 0{1/B). When 5^0, the maximum 
fluctuation of height of a solution F are of the order of 1/B. In this regime, 
we can show that with probability close to one, we can find interfaces that 
are suitably close to any smooth curve without self-intersections that connect 
(following its orientation) {0} x [0, L] to {L} x [0,L]. 

In order to understand the probabilistic problem we need to analyse the elemental 
components that an interface F can have. This is done in Section [2] where we analyse the 
solution curves of ([T]) and characterise them in terms of a free parameter that controls 
the height reached by the curve and its curvature. In fact, equation (jl]) implies that the 
curvature of an elemental component increases in absolute value with the distance to the 
reference height. The main results and their proofs are presented in Section [3l where the 
regimes outlined above are studied separately. Finally, Appendix deals with the proof 
of a technical lemma that is needed in Section 13.31 


2 Preliminary Analysis: Gluing of Solutions 


In this section we analyse the basic features of elemental components of an interface F. 
We first observe that equation ([I]) in the absence of boundary conditions can be solved 
explicitly in terms of elliptic integrals, see [13]. One can achieve this by first introducing 
the variable X 2 = X 2 — A and using the angle f3 G [— vr, tt] from ei to the tangent vector 
t = (x)^(s), 0 : 2 ( 5 )) where s is the arc-length. Then multiplying ([1]) by 0 : 2 ( 5 ) and integrating 
once we get 

—2 

fi^ + cos/3 = 0, (5) 

where O is a constant of integration that can have either sign or be zero, but clearly 
C > —1. Then solving ([T]), with given initial conditions, amounts to solving 

o:'^( 5 ) = cos/9 2 : 2 ( 5 ) = sin/3 (6) 


subject to the same initial conditions, together with (jSj) and a compatible constant C. 

The are the obvious symmetries of equation ([T]), ( 0 : 1 , 5 ) —?■ (2a —xi, — 5 ) (for any a G M) 
and X 2 —X 2 , that will be used to simplify the computations below. 

Using ([6|) and the definition of arc-length, one can write ([5|) as 


Xr 


2 \ 2 




2 


x: 


{dx,Y =[C-B:^] {dX 2 ) 


- \2 


(7) 


The fact that we are taking squares to go from ([T|) to o means that the set of solutions 
of the latter is larger than that of the former. In the construction one has to impose that 
the resulting curves are indeed solutions of ([I]), as we shall see. We also notice that all 
solutions must have 1x21 < ^2(U -|- l)/B. 
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Below, we find the solutions of ([1]) by gluing together pieces of curve solutions that are 
graphs of a function of the form Xi{x 2 ) and that solve ([7j). We then integrate ([7]) after 
taking the positive or negative square root of {dziYj{dx 2 Y] we use the equations: 


dxi 

dX2 



( 8 ) 


Observe that /3, for such graphs, is characterised by 

dxi 1 
dx 2 tan {3 

There are solutions, however, that cannot be constructed the way explained above. 
These are the solutions with X 2 = 0, that is the curves (s, 0) (air is on top of the liquid) 
and (— 5 , 0 ) (air is below the liquid). The former corresponds to a solution of ([5]) with 
0 = 1, and the latter to a solution of (jS]) with O = — 1. 

Next, we identify the points of the solution curves where either x'^(s) = 0 or X 2 (s) = 0, 
i.e. the values of X 2 where either the left or the right-hand side of ([7]) vanishes. For 
simplicity, we use the following notation; 

r /2(o-i) ^ ^ 

xT^{C-B)-.= l - if 0>1 (9) 

I 0 otherwise 


and 

ir'(C;B):=\P^F (10) 

Lemma 2.1 The points at which either x'^(s) = 0 or X 2 {s) = 0 are as follows: 

(i) If C > 0, when x[{s) = 0 then X 2 = B), and when X 2 {s) = 0 then either 

X 2 = ±a;™"'( 0 ;i?) or X 2 = ±a;^“^( 0 ; 5), except if0<C<l, that X 2 (s) = 0 only 
occurs at X 2 = ±a;™“^(0; B). 


(a) If C <0, then there is no value of s such that x'^(s) = 0, and as before, X 2 {s) = 0 
occurs at X 2 = ±x^°‘^{C; B). 

The next lemma describes the way solutions of ([I]) can be obtained by gluing solutions 
of (the branch equations) ([ 8 ]) at points with either vertical or horizontal tangent vector. 

Lemma 2.2 (Gluing rule) For C > —1, given the initial condition (a;i(0),0:2(0)) = 
{a,x'ff°‘^{C] B)) with a G M and (a;'^(0), 0:2(0)) = (—1,0) for (Q]), this specifies the branch 
equation of ^ that the corresponding solution (QJ) satisfies. If this solution reaches a point 
with o:'^(s) = 0 then it can only be continued further by solving the same branch of On 
the other hand, if the curve reaches a point with 0 : 2 ( 5 ) = 0 then it can only be continued 
further by solving the branch equation with opposite sign. 
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Proof. A solution stops being a graph of a function of the form xi{x 2 ) when, seen as a 
curve, X 2 (s) = 0. It is easy to see that, this never occurs at a point with X 2 = 0. At those 
points we can further continue the solution curve by gluing smoothly with another solution 
curve that solves (IH]) with the opposite sign (notice that these curves, as graphs, have the 
opposite orientation, and hence the change of sign in (JHl) implies that they have the same 
curvature at that point). 

On the other hand, at points with = 0, the solution is still a graph, and in 

principle, one could glue it to a solution that solves ([ 8 ]) with the opposite sign. However, if 
this happens at a point with X 2 7 ^ 0, then H would change sign, which contradicts ([I]). If, 
this happens at a point with X 2 = 0 , then H keeps its sign, but this also contradicts ([I]). 

■ 

We now construct pieces of curve solutions from which one can recover complete so¬ 
lutions of ([T]) (i.e. dehned for all s G M) by the symmetries (xi,s) —)■ ( 2 a — xi, —s), for 
any a G M, and X 2 —X 2 - The next lemma gives these pieces for every C > —1 by 
distinguishing the cases C> 1 , 1 >C >0 and 0 > C > — 1 . 

Lemma 2.3 (The pieces) Given the initial condition (a;i(0), X 2 ( 0 )) = (a, H), 

with C > —1 and a G M, and (a;((0), 0 : 2 ( 0 )) = (—1,0) for (QP, the elemental component of 
the solution curve is obtained by integrating with the minus sign, up to a X 2 = 0 or up 
to X 2 = B). Namely, 

(i) For C > I, the curve starts with dxi/dx 2 > 0 and reaches X 2 = B) for 

some value of s, for which a:((s) = 0 and 0 : 2 ( 5 ) = —1. The curve continues with 
dxi/dx 2 <0 up to a point where X 2 = x™"'(C; 5) > 0 for some other value of s, for 
which o:'^(s) = 1 and x^i^s) = 0 . 

(a) For C = 1, the curve starts with dxi/dx 2 > 0 and reaches X 2 = \J~^j for which 

x((s) = 0 and 0 : 2 ( 5 ) = —1. The curve continues with dxi/dx 2 < 0 and has an 
asymptote at X 2 = 0, that is with 0 : 1 ( 5 ) —)■ cxo and 0 : 2 ( 5 ) —)■ O'*" as 5 —)■ cxo. 

(Hi) For 0 < C < 1, the curve starts with dxi/dx 2 > 0 and reaches X 2 = B) 

at some value of s, for which o:'^(s) = 0 and 0 : 2 ( 5 ) = —1. The curve continues with 
dxi/dx 2 <0 up to a point where 0:2 = 0 af some other value of s, for which o:((s) > 0 
and 0 : 2 ( 5 ) < 0. 

(iv) For — 1 < C < 0, the curve starts with dxi/dx 2 > 0 and reaches X 2 = 0 at some value 
of s, for which o:']^(s) = 0 and 0 : 2 ( 5 ) = —1, if C = 0, and x[{s) < 0 and 0 : 2 ( 5 ) < 0, 
otherwise. 

Proof. We only prove (i), thus, let us assume that C > 1. We hrst observe that due to the 
fact that the maximum value that X 2 can take is the initial one, the direction of the initial 
tangent vector implies that away from, but close to, the initial point dxi/dx 2 > 0. Then the 
curve solution can be seen as a graph of a function of the form 0:1 ( 0 : 2 ) that solves ([H]) with 
the minus sign. Then the solution 0:1 ( 0 : 2 ) is convex as long as X 2 G (x™"'(C; B), xf)°'^{C] B)) 


(see (ID. that implies signed?xi / dx^ = sign(a: 2 )). In particular, this curve reaches a point 
with X 2 = x^‘^^{C-,B) (i.e. with (a;'^(s), X2('5)) = (0) “1) some value of s). In order to 
prove that the curve also reaches a point with X2 = (i.e. with {x[{s), X2{s)) = 

(1,0) for some value of s), we construct a curve ending up at such a point, and by a similar 
argument this curve also passes through a point with X 2 = a;™*^‘^(C'; B). Then, adjusting the 
initial value of xi so that both curves intersect at such a point (which they do tangentially), 
and using the uniqueness provided by Lemma 12.21 imply that they are the same solution 
curve. 

All the other cases are straight forward by convexity of the solution curve. ■ 

We show in hgures 2(a) (for C > 1), 2(b) (for 0 < C < 1), [3] (for C = 1), 2(c) (for 
(7 = 0) and 2(d) (for —1 < (7 < 0) the pieces described in Lemma [2.31 The next lemma 
describes how complete curves are obtained from these pieces. 


Lemma 2.4 (Complete curve solutions) For each C > —1, the pieces of solution 
curves obtained in Lemma \2.3\ give rise to complete curves, as follows: 


(i) There are three additional pieces of curve solutions obtained by the symmetries Xi —)■ 
—Xi and X 2 —t —X 2 - 

(a) In all cases above the curves obtained can be glued together smoothly, translating as 
necessary in xi, at the points where ^2(5) = 0 , and at points with X 2 = 0 (choosing 
pieces with the same tangent vector at the point). 


Remark 2.5 (i) The complete curves that we obtain are open and periodic except for a 

particular value of C for which the complete curves are closed: by continuity, there 
exists a value of C E (0,1), say Cq (actually Cq ~ 0.6522, see such that 

the corresponding elemental component satisfies that a:i(0) = Xi{x'f)°‘^{C]B)). In 
particular, if C > Cq, the curve obtained in Lemma\2f^has, globally, increasing xi 
as s ^ 00 , whereas the corresponding curve if C < Cq has, globally, decreasing xi as 
s —)■ 00 . 


Physically, one can interpret the solution curves with C < Cq as describing a situation 
in which air lies below the liguid, and vice versa for the ones with C > Cq. 


(a) It is also worth mentioning that the complete curve obtained for C G (—1,0] that 
oscillates around X 2 = 0 does not have self-intersections and can in principle connect 
two grains very far apart (in the horizontal direction). 


The next lemma gives the range of (3 that complete curve have depending on (7: 


Lemma 2.6 (Range of (I) For the complete curves obtained in Lemma 2.f the values of 
the angle (3 range in different intervals depending on C: 


(i) If C > 1, (3 takes all values in [—vr, vr]. 
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Figure 2: This picture shows solution curves of ([5]) with C > 0 and C ^ 1. The dashed- 
dotted line represents the line X 2 = 0, and the dashed horizontal ones are the lines X 2 = 
B). The thick solid lines represent the pieces of curves as specihed in Lemma [T3l 
The remaining lines give solutions obtained by the symmetries X 2 —X 2 and xi —?■ —xi, 
see Lemma 12.41 Any other solution is a translation of the ones shown. Complete solutions 
curves to CDi are achieved by the smooth gluing of these ones at the points with 0 : 2 ( 5 ) = 0 


on the dashed lines, (a) shows the case C > 1 and (b) shows the case in which 0 < C < 1. 


According to the orientation shown in these cases the air lies ‘on top’ of the graphs shown, 
(c) shows the case C = 01 and (d) shows the case in which —1 < C < 0. In these cases, 


the air lies below the graphs. 


(ii) If C = 1, (3 takes all values in [—vr, 0) U (0,7r], unless the solution is the horizontal 
line, for which = 0. 

(Hi) If —1 < C < 1, (3 takes all values in [—tt, — 77 ] U [ri,7i] for some 7] G ( 0 , 7 r], and 7] 
increases when C decreases (in fact, rj = 71 if and only if C = —1, i.e. X 2 = 0 with n 
pointing in the direction of gravity, and rj = 7 r /2 if and only if C = 0). 

Remark 2.7 (i) We observe that by the determination of f3 that we have chosen, its 

value jumps from vr to —vr at gluing points with X 2 = x't^°‘^{C] B), and, similarly, its 
value jumps from —tt to tt at gluing points with X 2 = —xlf'°'^{C] B). 

(ii) We observe that given the contact angle a and the radius of grains R, the contact 
point at a grain is given by the angle pc measured from the horizontal placed at the 
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center of the grain to the contact point, and satisfies 


pc = a - fic 


( 11 ) 


where fic is fi at the contact point. 

Then Lemma \2.6{ gives the possible range of Pc for a given a In particular, if C > 1 
Pc covers one period and unigue solution exists between two given grains that are 
sufficiently close by adjusting C. If — 1 < C < 1, then pc G [a — 7i,a—r]]U[a+ri,a+7i]. 
This means that there are some restrictions on the contact surface for grains 

that lie near the horizontal if air lies below the interface. 

We give the solutions of (|5]) with C = 1 explicitly in the next lemma, see also Figure [3] 



Figure 3: This picture shows four solution branches of ([5]) with C = 1. One if then is, in 
this case, X 2 = 0 (dashed-dotted line). The complete solution curves to ([I]) are achieved 
by smoothly gluing at X 2 = or at X 2 = —5) (dashed line), depending 

on the sign, and translating in the direction of Xi copies of the ones shown. 

Lemma 2.8 The solutions of ^ with 0 = 1 are explicitly given by the horizontal line 
X 2 = 0 oriented with n = 62 , and by two one-parameter families of solution branches: 

xf (x 2 -, a) = -^F{y/Bx 2 ) + a, X 2 e ^ > 

,0^ u 

where the parameter a is a real number and F is given by 

F(Z) = i log Z 6 (-2,0) U (0,2). (12) 

Moreover, F satisfies the following properties: hmz-)— 2 + F{Z) = lim 2 ^ 2 - F{.Z) = 
0 and Yimz^Q± F{Z) = —oo. Observe also that F'{Z) = 0 at Z = ±\/2 and that 
\\m.z^- 2 + F'{.Z) = oo, \\mz^ 2 - F'{Z) = —cxo and\\mz^Q± F\Z) = ±cxo. Regarding con¬ 
vexity; F'\Z) < 0 with Z G (—2, 0) U (0, 2). 



x^{x 2 ,a) =—j=F{\fBx 2 ) F a, X 2 & 
V B 


Vb 
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Remark 2.9 (i) We notice that the functions xf and x-^ can be extended to X 2 = 

±2/'/B, with xf{±2/y/B;a) = a and x^{±2/'/B;a) = a, but the derivative at 
those points is not finite. We also observe that gluing together the branch xf and xf 
with X 2 e [—2/\/R, 0), and the same value of a we obtain a complete planar curve. 
Similarly, xf and xf with X 2 G (0, 2/\/R] for the same a give another complete curve, 
that can be simply recovered from the previous one by the transformation X 2 —X 2 ■ 

(a) We observe that the limit C* —?• 1 can be studied independently of the value of B. One 
can easily show that solutions with positive C, and either C > 1 or C < 1 converge 
to the solutions found in Lemma \2.8\ on compact sets of X 2 ■ 

The rest of the section is devoted to deriving estimates on the maximal horizontal 
length of a solution curve depending on the value of C. First we derive the horizontal 
length at points where the height is minimum or 0 to points where the height is or 

^max. 

^2 

Lemma 2.10 Given a elemental component solution as in Lemma \2.31 ' 

(i) //—1 < C < 0, then 

(a) If < C < 1, then 

■^h,(C) < xi(0) - X, (xr^Ci B)) < ift2(C). 

where hfiC) ^ 0 as C —)■ 0 and hfiC) = 0(— log(-s/l — C)) as C ^ 1~, for both 
i = 1,2. 

(Hi) If C > 1, then there exist non-negative functions hfiC) and hfiC), continuous in C, 
such that, 

{xT^iC; B)) - xi {xr^C; B)) < , 

where hfiC) = 0(— log \/C — 1)) as C ^ I'*', for both i = 3,4 and hfiC) = 0{l/y/C) 
as C ^ oo, for both i = 3,4. 

Proof. In order to prove the estimates we use Lemma [2.31 and integrate ([8]) over X 2 in the 
appropriate interval where a:i(a; 2 ) is monotone. Then we estimate the integrand (that is the 
right-hand side of ([8])) with the minus sign by using straightforward bounds on the numer¬ 
ator and on the denominator. Namely, we estimate the parabola C — B^ from above and 

below. We also write (^1 - (c - = (l - ' (l + ^ 
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Then we either estimate the first or the second factor using the bounds on the parabola. 
The resulting integrals are elementary. We outline the proof of (i)^ the rest is left to the 
reader. 

(i) Let us assume first that — 1 < C < 0. Then 




..jmm 


C -B- 




1 

. 2 \ 2 


dZ, 


therefore 




V 


2(C+1) 


dZ 


TT 


B 


\/l - C” Jo (^1 + _ B^) " v^25(l - C) ’ 


and this shows the estimate on the right-hand side. 
On the other hand, 


- ^ >- 

l-C j-'d^ dZ _ 1-C / y2(C+l)+2 t ^ ^ 

ddcTDJdW (l-C + Bf)’ ~ /B(C+T) 


One can check that the right-hand tends to (-s/2—1) /\fB as O —>■ —1, this implies continuity 
of this term for C G [—1, 0]. Its minimum there is achieved at C* = 0, computing its value 
gives the result. ■ 

We now give the estimate on the horizontal length from points at height to points 
at height x^°‘^. 


Lemma 2.11 For any C > 0 the following holds 

-JgfdC) < X, (xr*(C;B)) -Xi {xr^{C;B)) < 

where h^iC) and he{C) decrease for C G [l,-|-cxo), h^{C), ho{C) —)■ 0 as O ^ -|-cxo and 
h^{C) tends to a positive constant as C ^ 0. 

We now adapt these estimates to actual curves that solve ([I]) subject to a Young 
condition on the grains connected. Before that we distinguish two types of solution curves. 
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Definition 2.12 Let pc and Qc G he the first and second, following orientation, contact 
points of a solution curve 7. We say that 7 goes from left to write if p^'^ < and write 
'jir- Similarly, we say that 7 goes from right to left if pc > qc^\ and we write 7 ^- 

We also denote by f3{pc) the angle that t forms with the horizontal at the contact point 
PcEMf. With this we define the parameter 


d{vi,R,C, PcQc) 


2R + 


cos ft (pc) - cos ft {qc)\ 
{2B{C + 1))^ 


> 0 . 


(13) 


Moreover, we shall indicate with a super-index + that 7 lies above X 2 = ViL (since it 
then has positive curvature) and with a super-index — that 7 lies below X 2 = viL (negative 
curvature). When 7 crosses X 2 = viL (its curvature changes sign across X 2 = viL) we 
indicate it with the super-index c. 


Lemma 2.13 Let and ( be centers of two grains connected by a solution 7 of m with 
a Young condition at contact points pc and qc (following orientation) and that solves ^ 
with C > 1, then: 

(^)Ifl = llr, 

d{vi,R,C,Pc,qc) < ll^-CII <2/2 + 2 (14) 

(^l)Ifl = lrl, 

d{vi,R,C,Pc,qc) < ll^-CII <2/2+ • (15) 

Proof. We write D := ||,^ — Cll (euclidean norm), then clearly D = 2/2 + \\Pc — gdl and we 
can write: 

WPc - qcf = \qi^'^ - + \pf^ - 

and 

\\Pc-qc\\ > \pf^ 

Using (jS]) then 

B-^ -^— = cos ft (qc) - cos ft (pc) 

which gives the right-hand side of flT^ and ([15]). 

On the other hand if 7 = 7 ^^ then 

0 < Ip7-+I < 2\x,(xr'‘(C;B))-xAxr’'{C;B))\ 


and 


\pI^'> - qf^\ < <“^(0; B) - B) 


\/2 2 
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If 7 = 7rZ then 


0 < \p^c^-qi^'>\ <2\xi{x^^^{C;B))-Xi{x^^'^{C;B))\ 


and 






< 


V2 2 

y^y/CTT+y/C^ 


The rest of the proof follows from lemmas [2.101 and 12.111 and the fact that — log |1 — Cj 
1/^|1 - Cl as C ^ 1. ■ 

The following lemmas are proved in a similar way than Lemma 12.131 


Lemma 2.14 Let ^ and ( be centers of two grains connected by a solution 7 of m with 
a Young condition at contact points Pc and Qc (following orientation) and that solves ([^ 
0 < C < 1.- 


(%) If ^ = 'jir, then 


d{vi,R,C,Pc,qc) < Ib-CII < -TTVT- 

{B{1 - C)}2 


(a) If j = ■jri, then 


d{vi,R,C,Pc,qc) < lb - Cll <2R + 


71^ + 2 

B{C + 1) 


1 

2 


(16) 


(17) 


We recall that for — 1 < C < 0 the complete curve does not intersect, so the estimates 
in Lemma 12.41 have to be summed up depending on the number of pieces of the curve (see 
Remark 12.5L ii)): 

Lemma 2.15 Let ^ and ( be centers of two grains connected by a solution 7 0 / (ll]j with 
a Young condition and that solves ^ with —1 < C < 0, then it is necessarily of the form 
■jri and satisfies 


1 

d{vi,R,C,Pc,qc) < lb - Cll <2R + n 
where n is the number of elemental components of 'yri- 


^ + 4)y 


(18) 


3 The Probabilistic Analysis 

In this section we shall give the main results and prove some of them. Before we do this, 
we complete the problem setting by describing its probabilistic features in detail. 
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We observe first that the probability of having at least m centers in 1/ C [0, L]^, PiV), 
is computed by 

We shall denote by VLy{L) the set of all possible configurations of centers in [0, L]^ 
that are distributed homogeneously and independently. In particular, there are centers 
in [0,L]^, and we can make the identification Vt^[L) = ([0,L]^)^'^ . We assume that the 
grains have the same size, where their radius R <1/2. Notice that we allow situations of 
overlapping grains. Then, we assign the following probability measure to VLy{L)\ 

= jj^Wdik, with d^ = Y[d^k- ( 20 ) 

k=l k=l 

The probability space is A^,, ^i^), where A^, is the a-algebra of Borel sets of flu{L). 

For simplicity of notation, we introduce the parameter 

6 * = L max i?|- ( 21 ) 

that relates the three (non-dimensionalised) lengths. As mentioned in the introduction, 6 
is relevant in the cases where B ^ 1. In the other cases we use to simplify notation; it will 
indicate that the solutions depend on B, R and L. 

For every to G ffi/(T) we denote by F(a;, A, 9) the interface solution of ([T]) and the Young 
condition. We recall that such an interface connects the domain xi = 0 to the domain 
Xi = L and that F(a;, A, 9) is the union of non-intersecting elemental components that are 
solutions of ([T]) and the Young condition between a pair of grains. 

In particular, each component of F(ci;, A, 9) solves (jS]) with a particular value of C, and 
we denote by C'(F) the maximum C for each such F. 

We index the elemental components of any F(a;, A, 6 *) following its orientation and let 
J = {l,...,m-|-l}be this set of indexes if the interface connects m grains. We shall write 
that a center G F(a;, A, 9) with i = 1,..., m if it is connected (by the Young conditions) 
to the elemental components 7 i_i and 7 * of F(a;, A, 9). We also denote by Vi = ||,^j — 

We shall need a number of results. The first one is a consequence of the lemmas 12.131 
12.141 and 12.151 and will be used in sections 13.11 and 13.21 


Corollary 3.1 For every L, to E and Vq, there exists a constant Dq, such that if a 

T{to,X,9) satisfies that 

D 

71 

then there exists a constant K > 0 such that every pair of centers of grains joined by F, 
fi and fi, and contained in || 0:2 — ViL\ > satisfy ||.^i — fi\\ < -|- 2R. In fact, 

Dooc 


I ^ 0 /or T) > Do, 


F(a;, A, 0) n < I a;2 — viL\ > 
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In many of the probability results we will use Stirling estimates (see e.g. 0 )- Namely, 


n\ = y/27rn ( — 
. e 


with 


1 


12(n + 1) 

Thus we can write for any n, m eN with 1 < m < n: 


< r„, < 


1 


12n 


n 


{n — m)m 


nr(m)-^(n 


m 


n+m 


( 22 ) 


(23) 


yiir 

We have the following lemma: 

Lemma 3.2 Given uq, n, m G N, with no < m < n — 1 then 

f ^ 1 with x = —, W{n,m) = - log(n) 

\m J n \n J 2 \n J 2 

with 

'h(x) = —xlog((r) — (1 — x) log(l — x) , (p(x) = log((r) + log(l — x ), 

and C{n) < for some c > 0. Moreover, at m = n/2, W reaches its maximum 

W{n,n/2) = —(n+1) log(l/2) —(1/2) log(?7,) >0, attaining its minimum at eitherm = n—l 
or m = no (depending on the value of no and n) for n large enough. 

Proof. The estimates follow from fl23l) . First, let K{n, m) = exp(r„ — Tm — fn-m) < C{n) 
with 


C{n) = exp(: 


< exp(: 


1) < exp(- 


11 


12n 12m + 1 12(n —m) + l' ~ ■^'12n 12n + 1 ' ' ~ \2n' 

where we use the estimates on (l22ll . The estimate on the other factor of (|23ll follows 
by a calculus exercise: It is easy to check that W{x) = n\l/(x) — has a maxi¬ 

mum at a; = 1/2 if n > 3 with lF(l/2) > 0. Also since lim„c^oW{x) = and lim^.^! = 
there must be at least two relative minima, one in (0,1/2) and one in (1/2,1). It is 
easy then to check that W''{x) > 0 in (0,1/2(1 — \/(n — l)/(n -|- 1)) and W'\x) < 0 in 
(1/2(1 — ^{n — l)/{n + 1), 1), where xi{n) = 1/2(1 — \/{n — l)/(n + 1) and a:)(n)l/2(l -|- 
{n — 1) / {n + 1) are the inflexion points. This means that there are two minima one is 
attained in (0, xiiji)) and the other in (xr(n), 1). Since Xriji) > {n — l)/n for all n > 1 and 
xiiji) < no/n with 1 > riQ < n, the statement follows. ■ 

3.1 The regime 1 < L <C min{\/5, R~^} 

The subcases of this regime can be treated similarly. Namely, the following theorem holds. 

Theorem 3.3 For all e > 0 there exists 9^ such that for all 9 < 9^ there exists fig C 
flu{L) with jJiyiRte) > 1 — o,nd such that for any oj G fie, the solutions T{u,X,9) satisfy 
r(ci;. A, 9) = {x 2 — ViL = h{u)} with 

\h(uj)\ < max 2R^ + 2KvoLR^ , 

where the positive constant K is as in Corollarv \3.1\ 
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Proof. For simplicity of notation we introduce the parameter 

5{9) = max 1^, 2i?| + 2Kv^LR^ , 
that satishes 5(6') —)■ 0 as 6^ —)■ 0. 

We hrst compute the probability of having a conhguration uj without centers in the 
strip \x 2 — vqL\ < 6. We call the set of such conhgurations flo, then, clearly, 

( 2^\ i^L'^ 

1 — j —y 1 as 9 —y 0. 

On the other hand if there exist cu G Oq and a r(a;,A,6') such that all its points satisfy 
\x 2 — VoL\ > 6 then, by Corollary 13.11 every two consecutive grains, and joined by F, 
satisfy ||^j — < M{9) where M{9) = j^^^ + 2R (observe that M{9) ^ 0 as 0 —)■ 0). The 

number of grains contained in F is at least Nq = This implies that the probability of 
such an event tends to 0 as 0 —)■ 0. ■ 

3.2 The regime L = 0(min{\/5, R~^}) and L —)► oo 

We hrst observe that if we denote by Uq the set of conhgurations for which the ’horizontal 
solutions exist (no grains, therefore, are connected), then /iy(Fo) = 1 ~ hy(<^o) where cuq is 
the set of conhgurations with at least one grain in the strip around the horizontal reference 
height of width 2R, thus = 2R/L and fiu{Uo) = 1 — 2R/L. In the current limit this 

means that /iy(Co) ^ 1 as L —)■ oo. 

If we know let Ui be the set of conhgurations that have at least a solution F with exactly 
one grain connected to the walls, then n^iUi) < where Ui is the set of conhgurations 

that have at least one center at a distance smaller than or equal than 2/\fB + i? to the 
horizontal reference line (recall that the supremum of the distance to the horizontal that 
can be reached is precisely 2 /\fB for the solution with C = 1 diherent from the horizontal 
line). Then /iy(a;i) = 2{2/\fB + R)/L —)■ 0 as L —)■ cxo and thus also ^v{Ui) —)■ 0 as L —)■ cxo. 

In the next theorem we obtain that in general fiu{Uk) —)■ 0 as L ^ cxo for all /c > 0 and 
that an upper bound for such behaviour is given by a Poisson distribution. 

First let us make precise the dehnition of the sets of conhgurations f/^: 

Definition 3.4 Let Uk C LL^{L) he the set of configurations u for which there is at least 
one solution F(a;, A, 9) such that k = {fj E u : F(a;, A, 9) n 0}. 

Theorem 3.5 If k>l, there exists Nq < 2Dq9tj < vLf and L sufficiently large such that 

No , 

i=k 

Moreover, for all e > 0, there exists a > 0 such that for all L > L^. there exists 
fie C flu{L) such that > I — e, and that for all cu G fl^ the solutions r(u,X,9) 

satisfy |C(F) — 1| < e. 
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Proof. First, we observe that the probability of having at least k connections is the sum 
of the probabilities of having exactly i connections with i = k ... vL^. Then we distinguish 
between the configurations that have connections outside the strip around the horizontal 
reference line of width 2DqR and those that do not. We first show that the probability of 
having connections outside the given strip tends to 0. On the other hand, the probability 
of having an interface connecting exactly k grains that lies inside a strip of width 2DqR is 
less than that of having k grains in the strip of that same width. Then we compute this 
probability. We recall the notation N = z/L^ that we use throughout this proof. 

First, we treat the case R S> Let ujg C kly{L) be the set of configurations that 

have a solution F(ci;, A,6*), with oj G Us-, that intersects grains ^ with — viL\ > DqR. 
In that case, there are at least three centers E u <Z Us with ||^j — ^j|| < (2 + K)R 

and, by Corollary 13.11 — ^z|| < (2 + K)R for some K of order 1. The probability of such 
events is: 




N\ (L2)^-2 

3 } (L2)^ 


'Ui-i2\\<{2+k)R 


di2 


= ^(iV - l)(iV - 2)7r2(2 + KfR^ ~ 
0 


as 


'U2-m<{2+K)R 
L —)■ oo . 




Let ul be the complementary event of Ug- Then, ^ 1 as L —)■ oo, and this together 

with Remark 12.91 (ii) shows the second statement of the Theorem. 

We let Uk the set of configurations that have exactly k grains in a strip of width 2DqR. 
Then, the probability measure of Uk is 



Where we have used that, in the current case, 6 = RL. 

In a strip of width 2DqR and the maximum number of centers expected is proportional 
to 2i'LDqR = 2uDo6 =: Nq = 0(1). To see this, we let Up>No be the set of configurations 
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for which the number of particles in the strip of width 2DqR is larger or equal to Nq. Then 


j=No 


j 


2Do9 


~ir) 


Using Lemma [3.21 we can estimate each term (except the last one which is {2Dq6/ 
as follows: 


where C{i>LY is as in Lemma [3.21 and 


< C{i>LY exp 





W{vL^,x) 


- + 0 {x log(x) + (1 - x) log(l - x)) - ^ log(z/L2) 

+ vL?‘{x\og{2DQ9/LY + (1 — x) log(l — 2Dq9/L^) with 


X e 


[2^00 \ 

L2 ’ • 


Then, W{yL‘^,x) —)> —oo as L —)■ oo, in particular, W{yL‘^,x) ~ vL‘^x\og{2DQ9/L"^) as 
L —)■ oo if X = 0(1) and Wx) ~ log(i/L^) as L ^ oo if x ~ Summing up all 
the terms implies that fLu{Up>No) ^ 0 as L —>■ oo at an exponential rate. 

Finally, we write 


Uk — {Uk n Ws) u {Uk n — {Uk n Ug) u {Uk n n Up>nY u {Uk n n Up^p^Y 

where Up^pq^ denotes the set of configurations for which the number of particles in the 
strip is less than Nq. Then, 

+ f; ("f) (2^)’ (i - YtY ' ■ 

As we have seem the first two terms tend to 0 as L —)■ oo, and the third tends to 
^- 2 D,eu 1 {2D,9ur as L ^ oo. 

The case R <C is analogous but working on a strip of width 2Do/^/B. 


3.3 The regime L ^ \/B with B ^ oo and i? —> 0 

In this regime we shall prove that the probability of having configurations with solutions 
that separate from the horizontal reference line a distance of order one tends to zero in a 
regime where L is not too large compare to some increasing function of B. 

We prove the result in the case 1/\Yb <C R, the other cases can be shown similarly, as 
we remark at the end of this section. 

Let us give two auxiliary lemmas, that will be needed in the proof of the main result 
of this section. The first one is a technical lemma: 
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Lemma 3.6 Let b > 0 and a sequence {a^} such that 


1 b 

Ofc + 7:: 


< flfc+l < Ofe + 2: 


21 + Ofc 1 + Ofc ’ 

and Oo > 5 > 0. Then, there exist positive constants Ci and C 2 independent of k, such that 

C'2\/fc < Ofc < Ci^fk. 


The proof follows by induction, multiplying the inequalities by Uk and au+i and summing 
up the results, the constants Ci and C 2 satisfy Ci > 46(1 + 6 /C 2 ) and C 2 > max 3 .( 6 y^ + 
^/h^x + 6 )/ 2 (l + Cl^/x). 

The next lemma is a geometric result, which we prove in Appendix lAl 

Lemma 3.7 //r(a;. A, 6) is a solution with m components joining (m—1) grains, i^i,..., fm-i, 
ordered following the orientation of T{u), X,9), then 

/ K \ 

(24) 

for some positive constants Ki and K 2 independent of L, B and R. 

Next, we define the configurations that may allow solution interfaces that separate a 
distance of order one from the reference horizontal line. The centers of such configurations 
that might be joined by a solution interface must satisfy certain conditions, according to 
Corollary 13.11 We gather these in the following definition. 


Definition 3.8 We let the union of the horizontal reference line to the left boundary and 
to right boundary be denoted by 

Zl = { 3^2 = viL) U {xi = 0} , Zl = {x2 = viL) U {xi = L) . 

Definition 3.9 For given constants K > 0 and h > 0 of order one, and a given 6 in 
the current regime, we let Uh{K, Dq,9, L) C flu{L) be the set of configurations of parti¬ 
cles such that for every u G U{K, Dq, 9, L) there exists Nq < and an injective map 
: {1,..., A"o} —)■ {1,..., z/L^} so that the centers {'Cc/(n)}n=i...Aro satisfy the following con¬ 
ditions: 

(i) max„{|^®) -niL|} > h. 

(ii) 

1^(2) _ f(2) I < J_^_ 
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(ill) //max{|Qi) - viLl, |Q(;+i) -vi^l} > Do/^^B then 

1 K 

|Kg(n+l) —'SsHII < ~ 

^ ICw - ^1^1 


(iv) 


dist{^g<^i), z[) < L)oniin 

dist{^g{No), Zl) < Do min 


J _^_I 


(v) Moreover, they satisfy [2^ . 

We observe that condition (ii) is suggested by the statement of Corollary 13.11 and 
the values of in ([9]) and flTOj) . Condition (Hi) corresponds is also suggested 

by Corollary 13.11 Finally, condition (iv) implies that the sequence starts either near the 
horizontal or near Zf and ends either near the horizontal or near Zf. 

We can now prove the following theorem: 

Theorem 3.10 Let Qh C denote the set of all configurations uj such that there 

exists a r(ci;, A, 6) satisfying max — niL| > h for some h > 0 of order one. Then, if 
y/B L <C y/B^ogB as B ^ oo, —)■ 0 as i? —)■ cxo. 

Proof. Definition 13.91 implies that jjiyiDh) < Do,9, L)). Let us then get an 

estimate on n^iUh^K, Dq, 9, L)). 

We denote by := {(? : {1,..., m} ^ {1,..., TV} : injective} and for every m < N 
and every g G we let C Uh{K, Do,9, L) be the set configurations for which there 
exists a family of centers {^g(n)}n=i,...,m satisfying the conditions of Definition 13.91 Observe 
that these sets of configurations are not, in general, disjoint. Then 


N 

Uh{K,Do,9,L) = U U 

m=l geSm 


and we estimate its probability as follows 

N 


fi,{Uh{K,Do,9,L)) = 


1 


L2N 


n 


d^n 


N 


N 


. N 


- 

N 

m=l 


EE nM 


m=l gGSm \n=l 


JUS 


Lin 


Y.N(N-1) 


m=l 


Ar(iV - 1)... (iV - m + 1) 

L2m 


^ m p 

n 

Vn=l dUs 


{N — m + 1) 
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where we have reordered the first indexes by setting g{j) = j with j G {1 ,..., m} for every 
m, and in the last inequality we have simply counted the number of elements in each 
and also estimated the measure of the elements with n > m by the total measure L^. 
Now, we consider two types of configurations in each , the ones for which all points 
satisfy that 

1 ^^^^ — ViL\ > for all n e {1 ,..., m} , 

vB 


and the ones for which there exists at least one j with 
we have, due to the properties in Dehnition 13.91 that 



< In the first case. 


K 

ll^n-^n+i|| < —for ne{0,...,m}, 

VoV B 

where we abuse notation by referring to the vertical lines Xi = 0 and Xi = L hy .^o and 
respectively, and denoting the distance of points to these lines with the norm symbol. 
Observe that then m > LDoy/B/K. 

Let us denote by 


and 

then 

with 


Uf = Wn{u;e 0,(L) : - v,L\ > ^ Vn e {1,..., m}} 


Dn 


Ul = W\{u:^ n,{L) : - v,L\ > ^ Vn G {1,..., m}} , 


fi,{Uh{K,Do,e,L)) <h + l2 




m=Mi 


L2r 


\\l_din. J = l,2, 

n=l“' 0 ' 


(25) 

(26) 


where Mi > LD^y/B/K and M 2 will be determined bellow. 

We hrst estimate Ji. For B large enough so that /{{Dq)^B) < 1, we have that 


N 


h< iV(iV-l)...(iV-m + l) 

m>LDo y/B j K 


L2" 


■ TT 


m—1 


K 


DoVB 


2m—2 


N 


m>LDo\^B / K 




2 \ "i-i 


.{DoYB 

and this tends to zero provided that 


< M 


LDq \/b 

utiK^ 5 K 


{D,YB 


(^)Lv^ ^ ^2 ^ 


(27) 
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In the current limit, ^ B and B ^ oo, then for large enough 6 it is clear that 
^ (5)^ and so holds and Ji —)■ 0 as L —oo. 

In order to estimate / 2 , for a a; G t/| we dehne a subset of the points in u as follows. First 
observe that, by the definition of 17|, there is a jo ^ {1; • • • > such that niL|| < 

on the other hand there exist a ji G {1,..., m} that —'^1-^11 > 2Dq{R+ 1/'/B). Then 
we can dehne the hrst point of the subset by taking and (the last point 

of the subset). For each A; G {1,..., M 2 }, we choose the points as follows 

then by the properties of Dehnition 13.91 and this dehnition, the subset satishes: 





v,L\ > Idf 

-lil < IC’ 


K 1 

ViL\ + 2 —-——- 

-KiL|) 


and also there exist a such that — niL| > Do/{2^/B). 

kQ 

We can apply Lemma 13.61 with = y/B\^l‘^'^ — viL\ to obtain an estimate on M 2 , 
namely. 


so 


h<Ci 

M2 > Mq '.= 

We then estimate the second term in 


Vf 

x/s’ 


h 


B. 


Cl. 
as follows. 


N 

l2< ^ 

m>Mo 

N 

S Z 

m>Mo 


N{N - 1)... {N - m + 1) 

L2m 

N{N -1) ...{N -m + l) 


L2n. 

'Mo 





Here we have indexed the elements of t/| \ by with n G {1,..., m — Mq}. Also, 

in the last inequality, we have used that the elements of satisfy ^ 

2K/B\^ii^ — viL\ and Lemma lT6l We recall that Vn = ||^n — ^n+i||- 
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Then, using the inequality of arithmetic and geometric means and Lemma 13.71 we 
obtain 


h < 


y/W \C2 ^/bJ 


< L" 


r) \ Mq—I ^ 

2K ^ ^ 

7n>Mo 

Mo-1 N 


m — 1 


m \ 


2 (m-l) 


n=2 




2K \ 


VW- \C2y/B 


' rm T\/T 


m>Mo 


^ 'm^ + K2L 


m-l\ y/B 


2(m—1) 


We rewrite 


u 


m—;= + K 2 L 


2(m—1) 


^m-l \ y/B 

then, for all m > Mq and B large enough 


m 


2(m—1) 


m K 2 Ly/B \ 
m — 1 Ki m — 1 ) 


2(m—1) 




I 2 < L 


2 (Mo-l) 


3 1 


y/W Ky/Bj 


We recall that Mq oc B and using fl2^ we can conclude that I 2 tends to zero as L —)■ 00 
provided that largest exponential growing factor is controlled by the fastest exponentially 
decreasing one, that is if Ly/B <C B \ogB as L —)■ 00 , and this implies the result. 


Theorem 13.101 can be formulated in the case 1/y/B < i? by replacing y/B by R~^ when 
i? —)■ 0. The proof is analogous, except for the proof of Lemma 13.71 In this case, the 
conclusion is simpler to achieve, by realising that the longest part of an interface joining 
the least number of grains is one crossing X 2 = ViL, all other parts are away from this one 
a distance of order l/yfB, applying Corollary 13.11 the lemma follows. 


3.4 The regime L = 0{1/B) as B 0 

We only consider here the case L = 0{1/B), and to simplify we assume that L = a/B for 
some a > 0 of order one. The reason for this is that we expect the maximum difference in 
height in a full solution T to be of order \/B. Indeed, if T reaches distances from X 2 = X 
much larger than 1/B, the elemental components of T on that region would have a very 
small radius of curvature, and could not join grains at an average distance of order one. 

In this section we further restrict the height of the possible T’s to the interval [A — 
EqL, X + EqL] where, essentially 

0 < £0 < min jui, 1 “ ^ 

and (for reasons that become apparent below) ^ Q- This ensures that the curvature of 
T in this region is small enough and will allow to approximate the elemental components of 
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r by straight lines, since for X 2 in this interval \H\ < aeo <1/2. We recall, however, that 
this restriction depends on the conhguration of grains, since A = ViL depends on F. We 
shall see that we can £x a volume V that is close to vi for the conhgurations of interest. 

We show that, in this regime, given any regular curve without self-intersections con¬ 
tained in [0, L] X [X — EoL, X + EqL] there are configurations co G with probability close 

to one, such that the interface solution r(a;,A,6*) is close to such curve in some suitable 
sense. 

Before we make precise this theorem, let us give the main idea of the proof and establish 
some settings and dehnitions. We make a partition of the domain into small squares Q. For 
any given smooth curve A without self-intersections and compatible with the prescribed 
volume condition, we show that the probability that a solution interface F satishes FflQ 7 ^ 0 
for every Q that has Q fl A 7 ^ 0 tends to one as L —)■ 00 . This means that all curves 
compatible with the hxed volume condition can be approximated by an interface F. The 
proof of that the probability ofFn( 57 ^ 0 ifQnA 7^0 tends to zero, will be formulated 
as a percolation problem, see below for the details. 

We can now state the main theorem of this section, but first we define the notion of 
compatible curves. 

Definition 3.11 (Compatible curves) Let V G (0,1) and Eq G (0,1) fl Q with 

2 a j 

We say that a curve without self-intersections A : [0,1] [0,1]^ is compatible with 

{V,Eo) z/A([0,l]) C [0,1] X [l^-£o,C + £o], F(0) G {0} x [V - eo,V + e,], F(1) G 
{1} X \V — Eq,V -|- £ 0 ] o,nd the area under F is equal to V. 

We remark that we compare real interface solutions F to compatible curves in the 
rescaled domain [0,1] x [C — Eq, V -|- £o]- We shall denote by F the result of expressing F 
in the variables x' = xf L. Then: 

Theorem 3.12 Assume that B = a/L for some a > 0. Given V G (0,1) let Eq as in 
Definition \S.ll\ and let A G C^([0,1]) be a curve compatible with {V^Eo). Then, there exist 
1^0 > 0 and Lq > 0 such that for all n > vq and L > Lq there exists > 0, such that 

oimax{R) 7 r /2 as R ^ t), and U G Tt,j{L) and a e{L) > 0 with B << ^(T), e{L) 0 as 
L —)■ 00 such that 

> ^e{L) with 0 < —> 1 “ as L —)■ 00 , 

and that for any oj G Ul, there exists a solution interface T{u,X,9) with contact angle 
a G [—'K/2,amax{.R)] that satisfies 

f (cj. A, 6 ) C T^^(^)(A), A C (cj. A, 6 )), 

with vi ^ V as L —)■ 00 . 


£0 < min < y, 1 — y. 
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Here stands for the tube of radius \/2e around a curve. That is it contains the 
points of the plane that are at a distance smaller than or equal to \/2e of the given curve. 

Remark 3.13 We remark, that ifvo G (0,1) is a given volume fraction of liquid, then one 
could tune the theorem a bit to obtain solutions for this prescribed volume with probability 
close to one. This can be done by taking V = Vq/{1 — i/irR) and R sufficiently small, since 
the expected number of grains under A is oL'^V and Vq is related to Vi by 

Next, we dehne partitions of [0, L] x [{V — eo)L, (V + eo)L] into small squares of equal 
size. In each of such squares that intersect A we formulate a percolation-like problem. 


Definition 3.14 Let 

Q:=[0,L]x[{V-eo)L,{V + eo)L]. (28) 

Let e > 0, depending on L, such that R -C £ -C 1 and such that 1/e, ‘le^/e G N. 

We define {Q}q£I^ to be a partition of Q into squares of size eL x eL such that if Q, 
Q E Is then QP\Q = ih if Q ^ Q, unless Q and Q have one edge or one vertex in common. 
We let M := e{L) L for simplicity. 

Clearly, = ‘Ieq/e^ < oo. We also observe that this definition requires ^ Q- 
We now consider a subdivision of each Q in smaller sub-squares of size of order one. 
We assume that the size is exactly one, for that we either need to assume that L G N or 
to rescale the domain. We henceforth adopt the first choice without loss of generality. 


Definition 3.15 Assume that L G N Let e{L) and Q E Is as in Definition 3. If. then we 


let Q = U^ejS'K where is a partition of Q into closed squares of size one that cover 

Q and such that their interior sets are disjoint. 

Moreover, for a s E {0, 1) and for each n E J, we let R^ <Z be a closed squared with 
\RJ = s centered in W. 


Observe that increases as L —)■ oo, but IS'^I = 1, and this value equals the average 
distance between centers (and is of the same order of the minimal radius of curvature in 
this domain). Then we can guarantee that if there are grains in neighbouring S'^’s, then 
there is a connection by a solution 7 of ([T]). In this way the proof reduces to estimating 
the probability of finding grains in the S'^’s that can form a connection between the given 
portions of dQ that intersect A. This problem is reminiscent of a site percolation problem 
in two dimensions (see e.g. na, HI), although the probabilities assigned to the sites, as we 
explain below, are weakly correlated. 

In fact, we shall restrict the connections to grains that fall in adjacent RAs. This is 
because we have to guarantee that elemental components do not intersect each other, this 
is possible by reducing IR^I = s, as we shall see. But first, we have the following result that 
guarantees the existence of 7 ’s joining grains in adjacent sites. This result is independent 
of £. 
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Lemma 3.16 (i) Let and S^' such that S'^HS'k' 7 ^ 0 with one edge in common. Then 

if there exist E Rk and f G Rk' then, there exists a Lq > 0 such that for all L > Lq 
there exists a solution of (QJj that joins and 

(a) There exists sq G (0,1) and amax{R) with amax{R) —^ 7 r /2 as R ^ 0, such that if 
K, k' and K G {1, ..., have a vertex in common the arc that joins a grain with 
center in f E R^ to a grain with center in G R^,' do not intersect the arc thar joins 
the grain Br{^') to a grain with center f E R^. 


Proof. Proof of (i): Observe that the centers f and f satisfy that ||^—f|| < Vs + 2s^ + 1 < 
2. On the other hand, \H\ < aso, in particular the radius of curvature of a solution joining 
these grains, say r, has |r|oo > 1/aeo > 2 . 

If there is a solution 7 that joins dB{f^,R) and dB{^,R), then the variation of the 
curvature of 7, AH < a |^2 — ^ 2 !/-^ < 2a/L (where AH = \Hmax{'l) — Hmin{T)\ is the 
difference of the maximum curvature and the minimum curvature of 7). This means that 
such a 7 is very close to a circle with large radius of curvature. 

On the other hand, there are straight lines that join dB{^,R) and dB{^,R) with any 
prescribed contact angle. Indeed, there are two parallel lines that are tangential to dB{^, R) 
and dB{^,R). One has contact angle 7 r /2 and the other has contact angle — 7 r/ 2 . Then 
for any given contact angle, there is a segment parallel to these ones that joins dB{^,R) 
and dB{^, R) with that contact angle. This is also true for circumferences of large enough 
radius, being the smallest such circumference the one that has the two grains inscribed at 
the largest possible distance; i.e. one with diameter ||^ — ^|| + 2R. 

We recall that H{s) = (3'{s) where f3 G [—tt, tt] is the angle from ei to the tangent vector 
t, thus |/9'(s)| < ae < 1/2 and |/3''(s)| < a/L. This implies that for L sufficiently large we 
can hnd a solution 7 close to a large circumference that joins dBijf, R) and dB{^, R) with 
a prescribed contact angle by adjusting the value of C and the height X 2 (see ([5])). 

Proof of (a): Let us see that reducing if necessary the size of each the solutions 
have no self-intersections. We concentrate in three adjacent sites, say (1,1), (1,2) and 
( 2 , 2 ) for definiteness and assume that there is exactly one grain in each of R{i,i), R{i, 2 ) and 
R( 2 , 2 ), with centers ^( 1 , 2 ) and ^( 2 , 2 ), respectively. Other situations can be recovered 

by translation and by a 7 r/ 2 , vr or 37 r /2 rotation. 

Let 6 s = (1 — \fs)/2 and let (fs be the angle that the segments that join the grain 
in R(i, 2 ) to the other two grains form. The length of such segments is \/s -|- 1 G (1, \/2). 
Self-intersections may occur in several ways, we aim to avoid then in the situations where 
the angle between joining elemental components if narrowest. The minimum possible (fs-, 
if s is fixed, is achieved when the centers of the grains are (1 — 6 ^, 1 — 6 ^) for the one in 
(6s, 2 — 6s) for the one in i?(i, 2 ) and (l-|-6s, l-|-6s) for the one in R{ 2 , 2 ) (see Figurell]). 
This gives: 


cos 0s 


2(1 - 26s) 

(1 ~ 262)^ -|- 1 


2^/i 

s -M ‘ 


(29) 


This function is an increasing concave function for s E (0,1), this implies that 0s £ (0,7r/2) 
is a decreasing function with 0s = 7r/2ifs = O and 0s = 0 if s = 1. 






First we assume that i? = 0, in that case we do not have to worry about contact 
angles. Suppose now that the grains are joined by elemental components that are circular 
arcs tangent at {bs,2 — bs), i.e. at the grain in R{i, 2 )- This gives the minimum radius of 
curvature 

/ X _ Vs + 1 

f'min\S) ~ r, ■ ( 1 /r)\ ■ 

2sm{(ps/2) 



In view of fl2^ . we obtain that rmin{s) is an increasing function with rmm(O) = l/\/2 and 
lim^^i Tminis) = oo. This implies that if the radius of the circular arcs is smaller than 
they intersect, i.e. this implies that for each s G (0,1), if 


H\ < min{2sin(0s/2)/\/FTT, 


(31) 


the arcs do not intersect. In particular, if we take 0 < s < Sq where Sq is such that 
^"mm(so) = this cnsures that \H\ < e^a and non-self intersecting arcs. 

Suppose now that R > 0, with R < 1/2. We can assume that H is nearly constant 
for hxed 0 < s < 1, so that fl3T]) is satished. Let r denote the radius of the corresponding 
circumference and (ci,C 2 ) the coordinates of its center. 

The situation is symmetric with respect to the diagonal with negative slope of 5(1^2) • 
The arc that joins i?(.^(i^ 2 )) to 5(^2, 2 )) (following orientation) has it center on the line 
X 2 = ;^(xi — 1) -|- 2. Let us assume, for dehniteness, that H > 0 at the contact point; in 
that case, ci > 6s -|- 1/2. We recall that if pc < — 7r/4, then, by symmetry, the arc that 
joins i?(,^(i_ 2 )) to B{^( 2 , 2 )) intersects the arc that joins i?(,^(i_ 2 )) to i?(.^(i_i)). Thus we have 
to hnd a sufficient condition on s that guarantees that pc > —7r/4, this will also imply a 
restriction on a. 

We can derive the radius of curvature for which the contact point has pc = —n/i 
and « = 0. This gives 


(V 


V(1 + + (3 - s)^ - 4:y/sR 

4Vstan((/s/2) 


which is an increasing function of s. The value a = 0 gives the maximum of the contact 
angle for this radius of curvature. Then, if s G (0,So) with sq such that rmin(so) = (a£o)~^ 
we obtain that r > rmm(s), since \H\ < aeo- Moreover, for such r given, there is a 
restriction on a, namely a < amax(s, r, R) < 0 where amax{s,r, R) is the contact angle 
formed by the arc that intersects i?(.^(i^ 2 )) at the point with p^ = —7r/4. 


In order to formulate the problem we identify each Q with a set of coordinates in 
namely: 


Definition 3.17 Let e and M = eL as in Definition 3.14 Adopting percolation terminol¬ 
ogy, we refer to each sub-sguare and, without loss of generality, to of the partition 
defined in Deftnition \3.15\ as a site. 

We identify every Q G via a suitable homeomorphism, with the sguare [0, M]^ C 
We also denote by dQ the union of sites that intersect the boundary ofQ. 
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Figure 4: Schematic picture showing how decreasing s may avoid self intersecting arcs. 
In this case the centers of the grains are placed at the corners of which makes the 
narrowest angle between arcs. 


Finally, the sites of Q = [0,M]^ are identified with the coordinate of their top-right 
corner; if we identify with n = {i,j), i,j G {1, 2, 3,..., M}, this denotes a site with 
top-right corner in {i,j). Without loss of generality, k = {i,j) will be identified with R^. 
We also identify the sides of dQ as follows: 

ai = {l,2,3,...,M}x{l}, a 2 = {l}x{l,2,3,...,M}, 
a 3 = {l,2,3,...,M} x{M}, a 4 = {M} X {1,2,3,..., M}. 

In the following we dehne what is meant by Q is ’connected’ or ’not connected’. This 
depends on whether the sites contain grains or not, and whether these are connected 
throughout Q. 

Definition 3.18 Let u E FtfiL), where we identify u = then: 

(i) We say that a site n E {1, 2, 3,..., is open if there exists at least one j E 
{1,..., uL"^} such that E R^- 

(a) Similarly, we say that a site k G {1, 2, 3,..., M}^ is closed z/ca fl /?« = 0. 

(Hi) The following random variable determines whether a site is open or closed: 

f 1, ^f > 0 

t{k-,uj)=1 (32) 

[ 0, if = 0, 

here xa denotes the characteristic function for some A C [0, M]^. 
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(a) The domain Q divided into squares Q 
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(b) A square Q divided into sub-squares S,^ 


Figure 5: This figure shows a schematic picture of the division of the domain Q that we 
use in the proof of Theorem I3.12[ In (a) we show a curve A contained in Q (thick solid 
line). The dashed line represents X 2 = X and Q appears divided into squares Q of size 
{sLy with B e 1. The shaded Q’s correspond the ones that intersect A. In (b) 


we 


show an arbitrary Q subdivided into squares of size one. Each contains the shaded 
region with = s G (0,1). Here a portion of dQ that intersects T is highlighted. 


The following definition sets when open and closed sites are connected, see e.g. [H] for 
similar definitions in two-dimensional site percolation. 

Definition 3.19 Let u G we say that two sites Kq, k G {1, 2, 3,..., are 

(i) adjacent, and we write ko71{u))k, if dist^Ko, k) = 1 and t{ko,u) = t{k,u) = 1. 

(a) connected if there exists a finite set of distinct sites, ki,..., Kj G {1, 2, 3,..., 

such that every two such consecutive sites are adjacent, and nQ7l{u)ni and Kj7l{u)n. 

A set {ko, ■ ■ ■, i^j, 1^} with this property is called a chain. 

(Hi) ^-adjacent, and we write koTZ*{uj)k, if dist^Ko, k) < \/2 and t{ko,u) = t{k,uj) = 0 . 

(iv) ^-connected if there exists a finite set of distinct sites, ki, ..., kj G {1, 2, 3,..., 

such that every two such consecutive sites are *-adjacent, and Ko71*{u)ki and KjTZ*{u)i<i. 
A set {kq, Ki,..., Kj, k} with this property is called a chain. 

The next definition gives the concept of regions being connected and *-connected. 

Definition 3.20 Let Ai C dk and A 2 C di for k 7 ^ I and k,l = 1,2, 3,4 he topologically 
connected regions of dQ. Then, for a given u G Qiy{L), we say that 

(i) Ai is connected to A2 if there exists a chain that intersects them. 
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(ii) Ai is ^-connected to A2 if there exists a *-chain that intersects them. 


We give a probability space to the configurations of sites induced by 


Definition 3.21 We endow VLm '■= {0, where Q = {1, 2, 3,..., with the structure 
of a probability space. We define the probability fiM as follows, for any A G Qm we set: 

P'MiA) = E : t{Q]0j) E A]^ . (33) 


For example, let = {k G {1,..., M}^ is closed} then 

/iM(C) = E VLu{,L) : = 0}) = “ 

and let (9^ = {^ G {1,..., M}^ is open}, then 

= p,v{{u: E : t{k]u:) = !})=!- fiM{Cn) ■ 



(34) 


We observe, that in contrast to the classical site percolation scenario, the probabilities 
of two different sites of being open are not mutually independent, because for finite L the 
number of grains is finite, and that influences the presence of finding more or less grains 
in different sites. 

On the other hand, since |i?K| = s for all k G {1,..., M}^ then fiMiCn) = q{L) G [0,1] 
and fiMiC^) = 1 — q{L) are independent of k G {1,..., M}^. Also, from (|34)) we get 




as L ^ 00 . 


(35) 


Thus, in the limit L ^ 00 , we find a classical site percolation scenario; in this limit 
the number of grains in the domain (as well as the domain size) goes to infinity, making 
independent of each other the events {u E flnu{L) : r(«); cn) = 1}. This means that there 
is a threshold value of u, say uq sufficiently large, such that the probability of finding an 
infinite cluster of closed sides is zero, i.e. that the dual event, finding an infinite cluster of 
open sites, has probability one, see e.g. H, HD], m- 
In our setting, we have the following percolation result: 


Theorem 3.22 (Percolation result) Let Ai C dk and A 2 C di for k 7 ^ I and k,l = 
1, 2 , 3,4 be disjoint and topologically connected regions of Q of size Ni and N 2 , respectively, 
with Ni, N 2 G (M/4, M/2). Let C(Ai, A2) C VLm be the set of configurations of sites for 
which Ai is connected to A 2 . Then, there exist z/q > 0 and Lq > 0 such that for all v > 1^0 
and L > Lq for e{L) as in Definition\3.T?\ such that 


fLM{C{A2,A2)) 


for some K > independent of L. 
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In order to prove this result we show that the probability of the dual or complementary 
events tends to zero. We identify the complementary events with the aid of the following 
topological lemma: 

Lemma 3.23 Let Ai C dk and A 2 C di for k 7 ^ I and k,l = 1,2, 3,4 be disjoint and 
topologically connected regions of dQ. Let A^ and A 2 he the topologically connected com¬ 
ponents of dQ\{Ai U A 2 ), then Ai and A 2 are not connected if and only if A^ and A 2 are 
^-connected. 

Results analogous to Lemma 13.231 are classical in percolation theory and are related to 
the Jordan curve theorem (in this case the exterior of Ai U A 2 in Q is ^-connected if and 
only if Ai and A 2 are not connected in Q). Notice that the lemma is independent of any 
probability choice. The proof can be adapted from the literature, see e.g. Bollobas [3]. 

We can now show Proposition 13.221 

Proof of Proposition 13.221 In view of Lemma 13.231 and taking into account that Ai 
and A 2 have length in (M/4, M/2), the crucial percolation result that we have to prove is 
that the probability of a *-chain to reach a distance of order M tends to zero. 

We hrst observe that given n G {1, 2, ..., M}^, let Chain(K,j) C LIm with j G N 
and j < M^, be the set of site conhgurations for which there exist a *-chain of the form 
{n, Ki,..., Hj-i}, then 


/iM(Chain(K, j)) < 9M 1 - 


L2 


uL^ 


(36) 


This is because the probability that a given set of j sites are closed is 

T \ 

i-L 

L2 


and there are at most 9^ ^-chains with that length (since there are at most 9 sites ^-adjacent 
to any given closed site). 

Let Jh G {1,...,M} the maximum horizontal distance from Ai to A 2 and let G 
{1,..., M} the maximal vertical distance from Ai to A 2 . Then J* = max{ J/^, satishes 
J* G (M/4,M]. 

Let C*{A\,Af) C LLm be the set of site conhgurations such that A\ and A 2 are *- 
connected, then for L sufficiently large 

/ ■ \ 

fLM{C*{Al, A^2)) <M • (37) 

Then, taking v large enough we obtain that log(9) — vs < —1/2, thus from fimi we get 

fiM{C*{A\,A2)) <Mr[M-\)e-^ . (38) 
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We write J* = IM for certain value I G (1/4,1], then 

This and the dehnition of M (Dehnition I3.15P concludes the proof. 

■ 

We can now show Theorem 13.121 

Proof of Theorem 13.121 Lemma 13.161 implies that we can determine the probability of 
a T intersecting Q through portions of dQ that intersect A by computing the probability 
of hnding grain centers in the sub-boxes This Lemma also give the restriction in the 
contact angle. 

We let uq and Lq and s{L) be as in Theorem 13.221 and assume that u > uq and L > Lq. 
Let Jl(A) C Id such that if i G Jl{^) then Qi C Ie(L) satishes that Q* U A 7 ^ 0. For every 
i G Tl(A), let also Si{Qi), S 2 {Qi) C dQi such that Sj{Qi) fl A 7 ^ 0 for j = 1, 2 and with 
length in (M/4, M/2) as in Theorem 13.221 
Then, by Theorem 13.221 

i.e. for all Ui C compatible with C{Si{Qi), S 2 {Qi)), > 1 — 1/(^(A)L)'^. On 

the other hand 

ieJz,(A) 

Observe that l/e{L) < < l/e^(L), thus taking e{L) = L~p with 0 < p < 1 we 

obtain the result with 

6,{L) = 1 - . 
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Appendix 


A 


Proof of Lemma 




We divide the proof in several steps. First we recall that we are assuming 
L S> min{\/R, with B ^ 00 and i? —?■ 0 . 


(39) 
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A rough estimate on the total length of a solution T that joins m centers is achieved 
by thinking of an interface having the longest possible length. We show that this worst- 
case scenario is given by an interface F having as many as possible folds elongating nearly 
horizontally from one vertical boundary to the other. It is clear that the longest elemental 
components of such interface will occur near X 2 = viL and that their union connects very 
few grains, since for such pieces C is close to 1 (see Lemma [2.101) . 

Taking into account that the volume is hxed and that the grains are uniformly dis¬ 
tributed (and do not see gravity) there must be folds on both sides of X 2 = viL. We can 
distinguish two types of interfaces, those that zigzag starting above X 2 = viL and those 
starting below it. 

Let us make the structure of such hypothetical interface solutions more precise. 

As we did in Dehnition l2.12l we also distinguish different types of elemental components 
and parts of the interface. First we distinguish those parts above X 2 = ViL, that we indicate 
with a subindex -|- (they have positive curvature), those below X 2 = viL, that we indicate 
with a superindex — since (negative curvature), and those that cross X 2 = viL, that have 
curvature changing sign across X 2 = ViL, and that we indicate with the superindex c. 
Then, we can dehne, skipping the dependency on (cu. A, 9) for simplicity, 

T = r+ u r- u T" with r+ = IJ 7+’* , r- = (J 7-’^ , and r" = (J 7"’'= , (40) 

i£j+ j&J~ 

and the set of indexes J’*', J~ and J‘^ are dehned in the obvious way, they are disjoint and 
J+ U J- U J" = J. 

We distinguish, also in Dehnition 12.121 two types of elemental components, those that, 
following its orientation, go from left to right (with subscript Ir), and those that go from 
right to left (subscript rl). 

We further decompose T into horizontal levels: 

Definition A.l (Horizontal level) Given an interface T, we say that a portion T C T 
is a horizontal level ofV if 

(^) 

min dist{fj, Xi = 0) < 2R + 1 and min dist{fj, xi = L) < 2R + 1, 


(a) And, if we let dminiX) = niiiixef dist{x,X2 = ViL) and dmaxiX) = miiia-gf dist{x,X2 = 
ViL), then 

dmaxX) - dyninX) < 2R + -j= —-^=--(41) 

V B 1 T yBdjjiiyiX) 

We now give a dehnition of a curve that zigzags horizontally: 

Definition A.2 (Zig-zag interface) A solution T of (QJ) with Young conditions is a 
zigzag interface if: 


35 










(i) There exists a set of consecutive levels Fq that contains F'^. One of these levels, that 
we denote by Fq, must ’cross’ X 2 = ViL, in the sense that if fc is the first and Cc the 
last contact point connected by Fg, then fc and (c on opposite sides of X 2 = viL. 

(a) There exists and N~ G N with N~, N~^ > 0 such that 



where F^ G F^ with i G {1,..., iV+} and F^- G F with i E {I,..., N } are horizon¬ 
tal levels ofV, indexed in decreasing order with proximity to X 2 = viL. 

Remark A.3 We distinguish four type of interfaces, namely: 

(i) Type I: the level Fq has elemental components of the form ■jir and crosses from positive 
to negative curvature. 

(a) Type II: the level Fq has elemental components of the form 7 ^; and crosses from 
positive to negative curvature. 

(Hi) Type III: the level Fq has elemental components of the form 'jir and crosses from 
negative to positive curvature. 

(iv) Type IV: the level Fq has elemental components of the form 'jri and crosses from 
negative to positive curvature. 

This definition does not exclude the possibility that other levels rather than the ones 
contained m Fq have elemental components of the form 'jir with C < 1. 

We shall write 


ro = 


U ) ^ ( U ^ ^ (U ) 


where, depending on the type of interface, Jq and Jq might be empty. 

Observe that there exists at most three consecutive levels in Fq; Fq, one above it Fq, 
and one below it, Fq^. The levels Fq and Fq ^ do not cross X 2 = viL in the sense that Fq 
does, but individual elemental components of them can. This is clear from the orientation 
of the elemental components and Dehnition I A. 21 (specify in at least one case of the types 
above). 

We might write: 





U 



for k = — 1 , 0 , 1 . 
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We shall first derive Lemma 13.71 for zigzag interfaces in which Lq has three levels. 
Then we prove that given the set of interfaces connecting m grains, then these interfaces 
maximise the length. 

Observe that such To’s connect at least 4 centers, a pathological case given by a Tq com¬ 
posed of just three elemental components, each being a level of Tq: Tq = {yj:;}, oscillating 
around X 2 = viL, Tq = just above X 2 = viL, and Tq ^ = { 7 ^ 7 }• 

Let us hrst derive some elementary properties of zigzag interfaces: 

Lemma A.4 If an interface T is a zigzag interface of type I or III and such that its Lq 
has exactly three levels, then: 

(i) and N~ are odd numbers. 

(ii) ^ 

rfi-i = U = U i = 1,..., ^ ^ ^ . 

(Hi) The index sets {J^} o.re disjoint with Ui=i...N+Jlt T and also {J)f} are disjoint 
with Ui=i...N-Jr T J~ ■ These sets might he empty. 

(iv) IfT connects m grains then 71 G and 7 m+i G ifV is of type I and 71 G Lj^. 
and 'jm+i e r++ if T is of type III. 

We leave the analogous lemma for type 11 and IV to the reader. Figure [ 6 ] shows a skecht 
of the four types of zigzag interfaces with three levels in Fq. 

Proof of Lemma 13.71 Assume there exists an ca G Tlu{L) and a F(a;, A,6*) that is a 
zig-zag interface and connects m grains where m > 3 and has three levels in Fq. it is clear 
any other interface that connects m grains is necessarily of the same length or shorter, 
the length of the interface in a strip of the size of Fq being at most of the order of 3L, 
outside the strip the length between connections of grains is of the same order as for a 
zigzag interface, according to Dehnition lA.ll (ii) and Corollary 13.11 

We can now prove the result for zigzag interfaces that join grains near X 2 = ViL with 
the maximal possible length. Let us assume that F is of type I. Then, it has three levels 
in Fq with three elemental components of the form 7 '^ whose length is of the order of L. 
Let 

no = if Jo , nf = ifJ+ and nf = ifJ- , 

where Jo = Jq H Jq" fl and jf are as in Lemma IA.4I Let us make abuse of notation 
by saying that ifJlj = 3, so that we index in this group only the pieces 7 '^ with length of 
order L, and the other pieces 7 ^^ of smaller order of length we index them in Jq is they are 
above the of the former, and in Jq if they are below. We then let = ifJ^. 

We recall that lemmas 12.101 and 12.111 imply that the other elemental components 7 '^ 
or 7 ^^ and in Fq have horizontal length comparable to the horizontal of elemental 
components in the even horizontal layers. 
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Figure 6: This picture shows a sketch of the four types of zig-zag solutions of Dehnition lA.2l 
The dashed line represents X 2 = uiT, the solid lines represent the interface. In the examples 
shown, the central horizontal level, Tq, has three parts (Tq, that crosses X 2 = uiT, Tq, above 
X 2 = uiT, and Tq \ below X 2 = V\L) each represented by a straight line. In the examples 
of type I and III there is just one horizontal level above and below Tq, and in the examples 
of type II and IV there are two above and two below. 


Thus we write 


P.+ P.+ 53 V,+ Y. P.' (42) 

2=0 2=0...m 2=0...m 2=0...m 2=0...m 

ieJQ iejy HJo 


Using Dehnition lA.il flTT|) . we can estimate the hrst term in fH2|) by 


Y, Vi <3 

2=0...m 
i&jg 




< 3 


-|- 2R -\- 



( 43 ) 


because dmini^ q) q ) 0. 

We write Vi = ((T>f)2 + (T>j')2)5 < Df+Dj' where V^ = and V^ = 
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In particular, 


1 = 0 ...m 


i£j, 


2=0...m 


i£j, 


0 


2=0...m 

i^Jo 


»r+ E ®."+ E D:<(N+ + N-+2)(2R+^y 


(44) 


Where we still need to control N~^ + N in terms of m. 

We let = minj=o,...,(jv±-i)/ 2 {’^^j} and = minj=i,,„,(jv±+i)/ 2 {nJ_J, and 

define = min{n^f”’+, n™"*’”} and = min{n;^f„+, It is clear that > 

2 and that rio > 4, then 

Af+ N- N+ N- {N+-l)/2 (Af--l)/2 

m+1 = n^+Y +^o = - (^^+l) + (^”+l)+ Y ^ 2 j+ Y ^2i+3 

2 = 1 2=1 2=0 2 = 0 j=0 j = 0 


and therefore, 


m 


+ l>iV+ + iV" + 5 + 


n 

I 


N+ + N- +2 


(45) 


Which in particular gives that m > N~^ + N +2. 

Let us get an estimate on the horizontal distances. Observe that 


E 


— '^min 


+ 1 


max V, . 


i=0,..,(7V±-l)/2 


ieJ, 


2 i=l,...,(V±-l)/2 

i^J2j 


2j 


(46) 


On the other hand, since the absolute value of the curvature of the pieces increases with 
the distance to the horizontal, for each j G {N^ — l)/2} we have X)iej±. — 

2E.E.,# ®?- Then, 


2j 


©.‘ss E 


is J2jUJ2j + l 


i^J2j 


We can now estimate the horizontal levels using (1471) and fl46l) . this gives 


(47) 


Y,Vi+ Y, ».'■+ E T’" < ^ max ©'■ + 2L, (48) 


ieJ^ ieJg 


is J^U J 2 J 


where the horizontal size of L N+ and of L n- is estimated by L. 
Using the estimates (143|) . (I44D . (I45|) and fl48|) in (|42|) we get 


^ / 4 \ 

Vi < (m — 2 ) ( 2R - — ] + 3(m — 4) max + 5L . 

^ " V /B/ OeJ 2 >J 2 - 


(49) 


39 









It is clear that the maximal distance between connected centers in the levels is 
achieved at an elemental component 7 with C = C* m an elemental component of the 
form 7^. Although we can only guarantee that C* > —1, we observe that the horizontal 
distances between centers are uniformly bounded. Namely, if C* > 1, Lemma 12.131 (ii) 
implies 


max Di < 2R + 
is 


( 71^ + 2 

[b{c* + i) 


<2R + 


f7r^ + 2 

V 2B 


1 

2 


If 0 < C* < 1 then Lemma [2.141 {(ii)) implies 


max Di < 2R + 

is J^U 


/ TT^ + 2 

V5(C'* + 1) 


<2R + 




1 


and if —1<C*< 0, then Lemma [2.151 implies 


max Di < 2R + 


|^ 2(7r^ + 4) ^ 2 


Thus, taking the maximum of these upper bounds we obtain that 


m 

^ A < (4m - 14) 

i=l 


( 2 R + 



+ 5L 


where K = \/2{'k‘^ + 4) 5.3 and this gives (1241) . 

The proof for the other types of zigzag interfaces is similar and is left to the reader. ■ 
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